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Abstract 

o . 

The process of optical frequency doubling can lead, in the undepleted regime, to the 

Q ' generation of a X-wave envelope with group velocity locked to the pump beam. Its 

^ . parameters and its angular spectrum, are directly related to the zero- and first-order 

C/3 . 

dispersive features of the nonlinear process. This constitutes a novel mechanism 
O i' for spatio-temporal localization of light. 
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The fact that the linear Maxwell equations admit localized solutions dates back to the 
book of Stratton, (l| as pointed out more recently. (0) However, this research topics has 
flourished only after Durnin et al. exjgeriment in 1987, who reported the practical realization 
of non-diffracting Bessel Jq beams. (3) Since then, linear waves localized in both space and 
time (i.e., non-diffracting and non-dispersing) have been reported in experiments, I^S;^), 
while main theoretical advances encompass the generalization of pro pag ation invariant solu- 



tions 



and the stud y o f their finite-aperture realizations [lol; lul; Il2)i or propagation 
in dispersive media (|l3t and waveg uid« U. The aspect of =uperluminality flH) a= 
well as potential applications such as femtosecond microscopy and communication systems 
have been also considered, (jiol : IitI i 

Recently, X-waves have been shown to be important also for nonlinear optics, since they 
exist as propagation invariant solutions to standard paraxial dispersive propagation models 
(3); whose formation is triggered by conical emission, (jiol ) In second harmonic (SH) gen- 
eration, these wavepackets are constituted by X-shaped envelopes at SH and fundamental 



frequency (FF) traveling locked together, which allow to explain recent experimental results 
()2Cl ). Since these X-waves continue to exist in the linear limit (in this sense they differ 
from any kind of soliton) where the FF and SH are uncoupled, it is natural to study the 
mechanism of their formation. In this Letter we show that a X-wave at SH can be gener- 
ated spontaneously via doubling of a broad FF pump beam in the undepleted stage of its 
propagation. 

To this end, let us start by reminding that the most simple axisymmetric localized wave 
in a dispersive medium with refractive index n{uj) and wavenumber k{uj) = ujn{uj) /c {c being 
vacuum light velocity) can be expressed as a superposition of Bessel beams (U) 

E{r^z,t)= / S{iu)Jo{Kr)exp[i{kzZ-~iut)]duj, (1) 
Jo 

where n = k{uj) sin[6'(u;)] and = k{u) cos[6'(u;)] are the transverse and longitudinal 
wavenumbers, respectively. (13) Equation(Q) describes a non-diffracting and non- dispersive 
beam provided that the following X-wave dispersion relation is satisfied 

k,{uj) =a + (3xi^, (2) 

Indeed, Eq.(j2I) entails vanishing dispersion dj^kz = at all orders n > 2, and in turn 
an intensity which depends only upon the retarded time t — Pxz, Px = K being the 



inverse wavepacket velocity. Given the spectrum S{uj), a and j3x are the two parameters that 
characterize the locahzed wave. In the case of vacuum, when a = the spatio-temporal 
spectrum has an X-shape {k = ±u;(c~^ — PxY^'^) and, due to invariant propagation, the 
far-field coincides with the near-field. Thus the wave retains its X-shape in the physical 
space {x,y,z,t) and turns out to be superluminal (c|/3x| < 1) (in the Fourier-transformable 
stationary regime no general theorem about causality and Maxwell equations is violated 



(21[). as discussed extensively in the literature (|22l: l23i l24l: 1251)). 



When a 7^ the X in the spectrum is curved and subluminal propagation is allowed. 



2fif ) In general these localized waves have infinite energy,^nd only approximations can be 
realized in practice (as for non-diffracting Bessel beams (3)). Solutions with finite energy. 



though no longer propagation- invariant, have been also reported. ()27l: |28|) In the case of 
dispersive media the situation is more complicated (the general angular spectra is given by 
= k{uy — {PxoJ + a)^, encompassing evanescent waves). However, strondv dispersive 
elements such as axicon and diffraction gratings or computer holograms (j^ permit to 
generate X-waves. 

Henceforth, we consider SH generation in the undepleted pump approximation. The 
paraxial and narrowband complex envelope amplitude a = a{r, t, z) at SH carrier frequency 
UJ2 = 2uJi obeys the following equation 

where AA; = k2 — 2ki {kj = k{ujj), j = 1,2) is the nonlinear phase mismatch, x the usual 
nonlinear coefficient such that |ap and are measured in W/m^, and e > stands for a 
loss term which is included for technical convenience, though we will take the limit e — > 0"*". 
We also neglect second-order dispersions k" = k"{ujj), since in practical cases their effects 
are usually small with respect to the dominant first-order dispersion measured by the group 
delay Pi — P2, Pj = k'iojj) being inverse group- velocities. The validity of Eq.Q requires also 
a broad (radially symmetric) pump beam so that its diffraction can be neglected, resulting 
into a pump amplitude A = A{t — Piz) = const., or in other words into a pump beam 
traveling undistorted at velocity /?f ^. 

Under these hypothesis, Eq. (jH)) can be easily integrated in the spectral domain by making 
use of the Fourier-Bessel representation: a{r,t,z) = a{K,Q, z)Jo{Kr)e~''^^KdK,dfl. 
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The solution of Eq.Q in terms of spectral variables (k, fl) reads 

. _ exp[i(^iO - Ak)z] - exp[i{(32^ - K^/2k2 + ie)z] 
~ - (32)^ -Ak + /cV2fc2 - le ■ ^ ^ 

where exp{iQ(3iz) is the Fourier-Bessel transform of the source term A^(r, t — f3iz) 

in Eq.(j2I), and we have considered, for sake of simplicity, the unseeded case, i.e. a{z = 0) = 

d{z = 0) = 0. Note that the denominator in Eq.(0]) represents (for e ^ 0) the expansion of 

the mismatch kz{2{uji + Q)) — 2kz{uJi + ^l). Consistently with our hypothesis, we limit the 

analysis to the case Pi ^ P2 (when Pi = (32, the effect of k" becomes relevant and will be 

addressed in the future). The inverse transform of Eq.Q in the lossless limit (e 0) can 

be obtained by applying residue theorems. We find a = whenever t < min{Piz, P2Z), or 

t > maxi^PiZ, P2Z), and 

t-p2z' 



a{r, t, z) = — — exp 



'2 1 



-iAk 



X 







00 . t-/3iz 

X{K)Jo{Kr)e^i^^KdK (5) 



elsewhere (see FigH for a pictorial representation of the support of a), with X = 
/2k2 + Ak)/{Pi — ^2),^^] determined by the pump spectrum. Hence in the region 
covered by the SH (travelling at speed 1/(32) which is emitted in the center of the FF beam, 
a is given by a non- diffracting beam which propagates at the group velocity of the FF. Indeed, 
Eq.(0) shows that the intensity |ap depends on z,t only through the combination t — (3iz. 
Only one-half of the localized wave is generated due to causality: the peak moves at the 
FF group velocity, while its costituents (the SH beam) moves at their characteristic velocity 
1/^2- As a result, the temporal extension of the wave is equal to \Pi — (32\z (see FigH]). This 
is similar to what has been theoretically predicted for practical finite-aperture realizations 
of vacuum localized waves. (Q) As in that case, it can been shown that the a field is the sum 
of two X-waves which cancel each other in the non-causal regions, a form of the extintion 
theorem. (3^ In the source-sink picture developed by Saari (17), only the source cone is 
present (as physically reasonable). 



Equation® strongly resembles the vacuum nondiffracting beams developed in (|26ri. How- 
ever, since we are dealing with a dispersive medium, we need to prove that the wave in 
Eq.(jni) belongs to the class given by (0). Consistently with the validity of Eq.®, we make 
the paraxial approximation kz = [k{ujY — k'^Y^'^ = k{uj) — /2k{uj), and then expand k(uj) 
around LJ2 at first-order, to end up (neglecting space-time couphng terms such as Qk^) with 



kz = k2 + (32^ — K^/2A;2, where Q = u — uj2. By equating this expression of kz to that in 
Eq.(j21), we can express the frequency as u = lj2 + {k2 — I3x0J2 — a — k'^ /'2k2)/ {Px — P2), which, 
after substitution in Eq.(IT)), yields 

t - hz 



E = 62 exp 



-i{k2 - [3x^2 - Ci] 



2 t-llx'^ 



/ X'{K)Ja{Kr)e^i^^^KdK (6) 
Jo 



where X\k) is a suitable redefinition of the spectrum Siu) and 62 = exp[i{k2Z — 002^)] is 
the optical carrier. The field in Eq.ijHj) has exactly the same form of that reported in Eq.©, 
once we account for the fact that the latter is an envelope which multiply the carrier 62- By 
comparing Eq. © and Eq. (0) we find that the parameters of the X-wave are directly related 
to those of the mixing process as 

px = Pi;a = k2-pxUJ2-^k = 2k^~P^LU2. (7) 

Therefore, our calculation shows that, for undepleted frequency doubling, the X-wave 
dispersion constraint Q is automatically satisfied and takes the form kz = PiQ + 2ki. 
Correspondingly the angular dispersion is given by the equation {k, = k29, 9 being the cone 
angle) 

which is an horizontal parabola in the angle- detuning plane, with position and convessity 
determined by the dispersive features of the process (group-delay and AA;). Its experimental 
observation would be a clear signature of the presence of angular dispersion and of the 
formation of SH-Xwaves. From Eq.® we can determine the energy of the SH S{z) = 
271 J°^^J'^\a\'^rdrdt = ^^^^^^£^z, where J'pF = J^\X\'^KdK. S grows linearly during 
propagation, its finiteness being a consequence of the X-wave compact support, and hence 
of causality. To make a practical example we assume the FF amplitude given by a gaussian 
pulsed beam with waist wp and duration tp- A = Aq exp[—{t — Piz^ /2tp — /2wp]. In this 
case J-'ff can be evaluated explicitly yielding 



fc2»pl/3i-/32l 



£ = ^-L^e- ^^'^ e-«[l + er/K)], (9) 
zy zntp 



with Sf the FF energy, ^ = k2wl{Ak - A(3x), A(3x = k2wl{(3i - (32Y/2t\, and a = 
tF/\/2\Pi — (32\k2VJp. Eq.© shows that for large |AA;| the energy conversion is negligibile, 



while the maximum efficiency is obtained for a particular value of A/c (in general Ak ^ 0) 
which can evaluated numerically. Choosing realistic parameters ()20l ) {wp = lO^^m, tp = 
lO^^^s, IPi — /?2| = 45 ps/m, k2 = 871 X lO^m^-*^), we have a = 1, and the maximum 
conversion efficiency is attained for Ak = A^x — — 1^21 /tpi i-e. the inverse spatial 
separation between FF and SH associated with tp. 

At last we wish to point out that different (often equivalent) terms are adopted to classify 
localized waves. Strictly speaking, Eq.(^ defines the so-called Focus Waves Mode (FWM), 
a more general class than X- waves, according to the literature. (j29|) However we are dealing 
with narrowband optical wavepackets, and hence FWM are known as X-Bessel (j^ Q). On 
the other hand Eq. (0) resembles the Focused X-waves (|l7|). Furthermore the causality 
implies that only one-half of the double cone is generated, and thus they can be identified 
with the so-called Y-, or Hankel-, waves. For the sake of simplicity we used the term 
"X-wave" from the beginning. 

In conclusion we have theoretically predicted that the process of optical frequency dou- 
bling can lead to an X-wave, traveling at the FF velocity, whose a parameter is directly 
linked to the wavevector mismatch. For a given input FF wave, the X-wave can be modi- 
fied by acting on the mixing control knobs (e.g. crystal temperature or orientation), while 
the resulting angular spectra (the far field) is given by horizontal parabolas [Eq.®]. This 
mechanism appears to be a valid alternative to the generation of optical localized-waves by 
means of linear techniques. 

This result is also susceptible of a different interpretation. The fact that the morpho- 
genesy into X-waves occurs spontaneously reveals the fundamental role that they can have 
in any ultrafast parametric process (e.g. four-wave mixing or optical parametric amplifi- 
cation) where the generalized phase-matching condition produces the correlation between 
spatial and temporal degrees of freedom which is behind the X-shape of the wavepacket. By 
studying such phenomena, the generation of new parametric sources could be conceived. 
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List of Figure Captions 

Fig. 1. Orientation (and support) of generated SH X-wave for Pi > P2 (top) and /3i < P2 
(bottom) . 
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FIG. 1 Orientation (and support) of generated SH X-wave for Pi > P2 (top) and Pi < P2 (bottom). 
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